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Abstract 

We employ the method used by Barbashov and collaborators in Quantum Field Theory to derive a path- 
integral representation of the T-matrix in nonrelativistic potential scattering which is free of functional 
integration over fictitious variables as was necessary before. The resulting expression serves as a starting 
point for a variational approximation applied to high-energy scattering from a Gaussian potential. Good 
K^ ' agreement with exact partial-wave calculations is found even at large scattering angles. A novel path-integral 

j^ ' representation of the scattering length is obtained in the low-energy limit. 



1. Path integral methods are not only useful in bound-state probleins but also give new insights and 
approximation schemes in quantum scattering. Recently two path-integral representations for the T-matrix in 
nonrclativistic potential scattering have been derived in Ref. [1] starting from the S'-matrix as (infinite) limit 
of the time evolution operator in the interaction picture. For a nonrclativistic particle of mass m moving in a 
potential V(x) they take the form 



m 



S{h) - 1 (1) 



where 



K == i(k,+k;) , i^=|K| = /ccos^, q = ky-k, , 9=|q| = 2fcsin^ (2) 

are the mean momentum and momentum transfer, respectively. E = fc^/(2m) is the scattering energy (we set 
h= \) and the scattering angle. 

The main features of these representations are functional integrations over velocities without boundary 
conditions and the use of "phantom" degrees of freedom to get rid of explicit phases which would diverge in the 
limit of large scattering times. Two versions exist which are distinguished by the reference path along which 
the particle dominantly travels and the dimensionality d of the "anti- velocity" w(t) which is needed to achieve 
the cancellation!^. 

S{h) =^ I V^vV^w cyi^li I dt^ [v2(i)-w2(f)] I e^xe^'^'^l (3) 

X(b, V, W] = / dtvih + Xrct {t) + Xquant (i, V, w] j (4) 



In the first case 



X 



{d=S.),,. _ K ^(d=3) 



Df 



(t) = — t, yi"- >.{t,w,^] = Xv(t)-Xw(0) (5) 



m 



the reference path is a straight-line path along the mean momentum and 

1 /■+°° 
x,(t) = -j di'sgn(t-t')v(t'), x,(i) = v(t) (6) 

describes the quantum fluctuations, sgn(a;) = 20 (.x) — 1 is the sign- function. In the second case the anti- velocity 
is only 1-dimensional 



■'^rcf \ I — ^ray(,tj 



^e(-i) + ^e(i) 
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(d=l) 



t^ X^uant(i,V,w] = Xv(t) -Xiv(O) -Ka:||„,(0) (7) 



and the reference path is a ray along the initial momentum for i < and along the final momentum for i > 0. 
In Ref. [T] it has been shown that both these representations reproduce the Born series to all orders if expanded 
in powers of the potential. 

Is it possible to get rid of the "anti- velocity" altogether? In the present Letter we show that this is feasible 
by using a method which has been introduced by Barbashov and coworkers [2l S |4] nearly 40 years ago. They 
showed how one can amputate exact Green functions in Quantum Field Theory without using perturbation 
theory but their work doesn't seem to have received much attention apart from the (rather specialized) field of 
quantum gravity where Fabbrichesi, Han and others [S] HI El [5] have extensively used the Barbashov method to 
describe gravitational scattering. Our main aim in this work is to transcribe their method to the nonrclativistic 



^Here and in the following the path integrals arc normalized such that the free (Gaussian) integral is unity. Our notation 
indicates that X i^ ^ function of the impact parameter b but a functional of v(t) and 'w{t). Similarly for Xquant- 



case and to use the result for an approximate variational calculation similar to the one investigated in Refs. 
[21 [To]. As a by-product a new path- integral representation of the scattering length is obtained. 



2. Our starting point is the following "reduction formula" for obtaining the on-shell T-matrix from the full 
Green function 



^f 



T 



k2=k2=2mB k2,k2.-i.2mB 



lim (kf (e-Ho) G(E) (e-Ho 



k, 



This is easily proved by using the definition of the exact Green function 

1 1 



GiE) := 



E - H + iO E- Ho-V + iO 



(8) 



(9) 



and some standard algebraic manipulations: 

E - Ho^ G{E) (^E - Ho^ = (^E - H + V^ G{E) (^E - Ho 



l + VG{E)\ iE-H + V 
= E-Ho + V [1 + G{E) V) = E-Ho + f. (10) 



When going to the energy shell the first term on the r.h.s. of Eq. ([TO)) vanishes and we obtain indeed the matrix 
element of the T-matrix between initial and final momentum states. Eq. ([5]) is the nonrclativistic counterpart 
of the standard Lchmann-Symanzik-Zimmermann (LSZ) reduction formula given in any textbook on Quantum 
Field Theory (sec, e.g. Eq. (7.42) in Ref. [H]). As other tools arc available, this relation seems not to be 
employed in standard nonrclativistic scattering theory. One exception is Ref. [12] where the time-dependent 
version in Eq. (16.82) is used to illustrate the "close analogy" of the reduction formula "with the S matrix 
constructed in propagator theory" . 

Assuming that the interaction V is given by a local potential ^(x) we now employ the Schwinger represen- 
tation of the full Green function © 



G{E) = -I /"dre<^-^+*)^ 

"'0 



(11) 



and the velocity path integral for the time evolution operator as given in Eq. (4) of Ref. [T] . Then one obtains 
from Eq. ^ % 



1^2 \ / k2 

%^f = i-i) lim [E -]\E ^ 

■* k2,k^^2mB V 2m} \ 2m 



dTe'^^ / (i3a;rf3^/g-»k^.x'+,:k..x 



X f V^v (5^3) ix -X- f dr v(t) j exp i i j dt —v^{t) -~V (x+ f dr v(t) 

Performing the x'-integration and shifting the velocity v(r) = v(t) -|- k//?n yields 
Ti^f = (-i) lim (A, A/)-/" dTe'^f'^ fd^xe-'^^f-^'^'' 

X / V^v exp <! i dt ^v'^{t) -V i x+ -^t + dr v(t) 



(12) 



(13) 



^Our scattering states are normalized to ( / | * ) = (^i')'' <5"' ( kj — k^ ) . 



where 



E 



2m 



^f 



2m 



(14) 



is the initial or final off-shcUness. It is obvious that Eq. (J13p is only non-zero if the Fourier-transformed Green 
function develops two poles when k^, k^ — > 2mE. How is that possible? It only can originate from the integral 
over T, presumably from T — >■ oo (one intuitively expects that the particle becomes on-(energy-)shell for infinite 
times). Indeed for any smooth function /(T) one has 



/•OO 

-i lim A / dTe^('^+*)^/(T) 



— lim 

A 



/"OO f) r 



(partial int.) 



/(O) + lim / dTe'^^ 



df{T) 
dT 



/(^). (15) 



However, integration over T only takes away the factor A/ but not the other factor A,; and one obviously 
needs a second integration. The main trick in Barbashov's treatment is a way to provide for this second time 
integration by first subtracting from Eq. p^ a term (a "1") which does not contribute in the on-shell limit 



%- 



-i lim A,Af / dTe''^'^f+'°^^ d^xe-'^i-'' 
/ V^v exp [ i / dt—w'^{t) j < cxp -z / d< F(p^(x,i, v]) 



X / P^w exp I i I dt^v''{t) I ■( cxp -i / d< F(p^(x,i, v]) - 1 )■ . (16) 

This is because the free Green function is proportional to J'^' (k/ — ki)/(_E — k^/(2m)) . We use the abbreviation 

Pj(x,t,v] :=x+-^t+ / dTv(T) (17) 

"^ Jo 



(18) 



and omit the bar on the velocity variable. Now the factor in curly brackets in Ec^. (|16p can be written as 



exp 



i / dtV{pj{-x.,t,v]) —1= ^* / rfC ^(P/(x, Cj v]) / dAcxp —iX / d<y(pj(x, t,v 
Jo Jo Jo Jo 



which provides an additional integration over the time ^. 

3. Since the algebra in Ref. [6] (which we follow with several corrections) is somewhat sketchy and involved 
we list in more detail the steps to be taken in the following: 



(i) Write 



dT / d^.. = dS, dT. 

Jo Jo J? 



(19) 



and substitute T =: ^ + s so that s £ [0, oo] . We then have 

%^f = i~if ^ hm^^A,A/ r d^ /""dse'(^^+*H«+^) f V^v exp [ z / dt'^^^t) j 

/(•I r ri+s 

d^xe-''i''V{pf{x,^,v]) / dA exp -iX dtV { Pf (x, i, v] ) 

Jo Jo 



ii) Substitute 



kf f^ 

X = Pf(x,^,v] = X+^^+ / dTw{T). 

m Jq 



(20) 



(21) 



Collecting the pieces in the exponents we then obtain 
Ti^f = {-if lim Aj A/ / d£, cxp 

Ai,A/-!-0 Jo 



*(A, + ^o)e + *g^e 



i+» 



X / ds cxp [i ( A/ + iO)s ] / P'^w cxp i / dr— v^(r)+iq- drwlr] 



/^i r /-C+s 

d^x e-*i* y(x) / dA exp -iX di F (x + py(x,i, v] - p^(x,^, v] 

■/o Jo 



(iii) The unwanted pieces (underlined in Eq. (j22p ) are eliminated by yet another shift 

v(t) = v(t) + -^ e(e - r) . 
The argument in the last potential term then becomes 

x + P/(x,t,v] -pj(x,^,v] = X- 

Finally, we define 



ke(t-e) + -e(e-t) 

m m 



(22) 



(23) 



(t-O+y drv(r). (24) 



f := T-e, i := i-^, v(f) := v(f + 
so that the argument p4)) becomes 



-e(t) + — e(-t) t+ / drv(f) =: p,,y(x,t,V] 
Omitting the bars over the transformed variables we thus have 



m 



(25) 



(26) 



%^f = (-i)' lim A, A. /d^e'(^'+'")« /" dse*(^/+^°)^ /d 

A.,A/^o ■* y Jo J 



^x e-'i-^ 



F(x) 



1 " 

X V^v cxp i dT^w'^ir) / dA exp ~i\ dtV {p^^y{x,t,v]] 

[ -^ J L -« 

It is now possible to truncate both external legs by using Eq. (|15p and to obtain the final result 
Ti^f = I d^xe-''i''V{x) jV^vexp i j dt^v^{t) j dA exp -iX j dt V {p,^y{x,t,v]) 



(27) 



(28) 



Of course, the A-integration can be performed trivially but for many applications it is better to leave it in this 
form (see below). 

Note that in the second potential term the particle propagates mainly along a ray formed by incident and 
outgoing momenta exactly as in the "ray" representation ([7]) while the quantum fluctuations are described by 
the functional integration over the velocity. Indeed 



/'rav(^'^'V) = X + Xray(i) +Xv(0 -Xv(0) . 



(29) 



Eq. (P5|) also gives a precise definition for a local "pseudopotential" (in general, this is a nonlocal object; see, 
e.g. chapter 6.2 in Ref. |13l ) whose Fourier transform yields the exact T-matrix. 



While it is not possible to set fc = directly in Eq. ([T]) the representation (|28p allows that without 
impunity. Therefore, defining the scattering length a as the negative (fc — ;• 0)-limit of the scattering amplitude 
/ = — TOT'/(27r) (as in Eq. (X.47) of Ref. [14]) we immediately obtain a path-integral representation of the 
scattering length 

a == ^ I (fxVin.) fv^vcxp if dt^v^it) I dXcxp -iX j dtv(x+ f drvirU . (30) 

Unfortunately, we do not have a simple argument why the scattering length should be real in potential scattering, 
i.e. why the imaginary part of the r.h.s of Eq. pop vanishes. This is one of the shortcomings of a path-integral 
representation of scattering; others are that basic properties like unitarity (and therefore the optical theorem) 
are more evident in the traditional operator formulation of Quantum Mechanics. 



4. We can use our result as a starting point for new approximation schemes. For example, following Refs. 
[9l [To] we may variationally approximate the velocity path integral in Eq. ((28|) 



F(x, A, S] 



V^v e*'^(''''*'-"^l 



S'(x, A, v] 



+ck: 



dt 



yv2(i)-Ay(p,.ay(x,i,v]) 



(31) 



by means of the Feynman- Jensen variational principle. The most general quadratic trial action, which allows 
the analytic evaluation of the needed path integrals, is 



SJv] = — v^Av + Bv 



(32) 



in the condensed notation of Ref. |10j where summation over discrete and continuous indices is implied. 
Evaluation of averages and derivation of the variational equations for the functions Aij(t,t'), Bi{t) follow the 
same lines as in Refs. [9j [10] where more details can be found. Here we just give the results: with this trial 
action the variational approximation to the T-matrix with no anti-velocity reads 



T 

' t- 



(AB30) 



«->/ 



d-'xe-'^^'^Vix) / dAi^var(x,A) 




i^var(x. A) 



JXo+iX^-Q. 



with 



^0 = 

(So),,(i,t') = 



dt U{t) 



-VUT^qVU 
2 



\t-t'\-\t'\-\t\ 



2m 
--Tr[ln(l + Ei7) -SiJ] 



Here U{t) is the Gaussian transform of the potential XV, i.e. 



U{t) 



and Hi 



(27r)3 



XV{p) exp ip ■ p^^,{t) - - p a{t) p 



with 



V{p) 



(33) 

(34a) 

(34b) 
(34c) 
(34d) 



d^x e-'P'' V{x) (34e) 



: didj U its Hessian. The width of the Gaussian transform is given by the diagonal value of the 
matrix S, i.e. crij(t) = Tiij(t,t) . Varying with respect to the variational functions A, B leads to variational 
equations for the trajectory and the (off-diagonal) width: 



PvarW = X + X,.ay(i) + (I]oV[/)(i) 



(35) 
(36) 



The first one is a classical equation of motion for the trajectory of the particle as can be seen by differentiating 
Eq. ((35|) twice with respect the the time t. This gives Newton's equation with a "kick" at f = 



mp^^M = q 



+ 00 



dt' VU{t') S{t) - VU{t) 



(37) 



The second one is a Lippmann-Schwinger-type equation for the quantum-mechanical spreading of the scattered 
wave. In the present numerical evaluation of these coupled nonlinear equations we have restricted ourselves to 
the case where A = 1. This amounts to setting the Hessian to zero and to neglect the wave spreading. Then 
the width of the Gaussian transform is fixed as Uij{t) = 6ij \t\/m and the quantity D, vanishes. 

This approximation (B30, in the nomenclature of Ref. [TU]) should reproduce the 2"*^ Born approximation 
and therefore be roughly as good as the systematic approximations studied in Ref. |10j which included the 
correction by the second cumulant. Indeed, simple algebra shows that the second Born approximation is 
contained in our variational approximation: When expanding T^_^f up to second order in the potential we 
can neglect the phase Xi and set exp(iXo) = 1 + iXq. Furthermore the argument of the potential in the phase 
Xq reduces to x + Xi.of(i) in this order. The A-integration is then trivial and with the definition of the Gaussian 
transformed potential one can perform the x-integration. The result is 



T, 



(B30) 



'^/ 



T 



1st Bor 



(27r)^ J_^ 



ip2-Xi.ay(t)-— |t|p2 



(38a) 

where the Gaussian width in the the "B30" -approximation has been inserted. This allows the t-integration to 
be performed and gives 



T, 



(3-0) 



i^f 



T 



lstBorn_i f fpiCPp_2 ,~, 

(27r)3 



V{pi)V{p2)S^^Hp,+P2~q) 



P2 • k//m - p|/(2m) + iO P2 • k^/m + p^/(2m) - tO 



0{V^) 



(38b) 



where the two terms come from the integration over positive and negative times, respectively. Now, putting 
Pi = p — ki , p2 = k/ — p in the first term and just reversed in the second, one sees that both terms are 
the same due to energy conservation E = k|/(2m) = k^/(2r7T.) which is already imposed by dealing with the 
on-shcU T-matrix. This cancels the factor 1/2 from the A-integration so that the result is 



rpiB30} _ rplstBorn 
^i^f - ^ 



''' Vi^f^p) 



1 



(27r) 



E - p'^ / {2m) + iO 



V{p~k,) + 0{V'') 



= T 



1st Bo 



+ k/ 



V- 



1 



E - Ha+iO 



V 



0{V'') 



3\ rrnlstBorn 



T 



■2nd Bo 



+ 0{V^) . (38c) 



We have solved the variational equations by iteration similar as was done in Refs. [HI HH] apart from 
numerical improvements to speed up the tabulation of the trajectory and to add the asymptotic contribution to 
integrals with slowly decreasing integrands. Fig. [Ijshows the comparison with an exact partial-wave calculation 
for scattering from a Gaussian potential 

V{x) = Foe""'/^' (39) 

at kR = 4, 2mVb i?^ = —4 which served as a benchmark calculation in Refs. [Hllin]. Plotted is the relative 
deviation |A///| = | ( /var ~ /exact ) //exact | and it is seen that our variational approximation outperforms 
both the first and the second Born approximation. Even at large scattering angles where the cross section is 
down by many orders of magnitude the deviation is still rather small. 
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Figure 1: (Color online) The relative deviation | A///| of the variational amplitude in the present approximation 
("B30") from the exact scattering amplitude as a function of the scattering angle 9. The error bars denote the 
estimated uncertainty of the numerical evaluation. Also shown arc the results from the first and from the first 
+ second Born approximation. 



5. In summary we have presented a new path- integral representation for the T-matrix in nonrelativistic 
potential scattering based on an ingenious method for amputating full Green functions in Quantum Field 
Theory due to Barbashov et al. Compared to previous representations the present one does not need artificial 
"phantom" degrees of freedom to cancel unphysical contributions. 

Of course, there is an (infinite) variety of path-integral representations depending at which level of formal 
scattering theory one introduces the path integral. For reasons of classification one may count the powers of 
the potential in front of the velocity path integral. Then the previous forms displayed in Eqs. ([1]), ^ may be 
characterized as "level 0"- representations whereas the Barbashov-like result ([25]) belongs to the "level 1" -class. 
A path-integral representation based on the formal solution of the Lippmann-Schwinger equation (last line in 
Eq. (fTO|) ) would lead to a "level 2" -representation and so on. Obviously, higher-level representations become 
more and more complicated so that the formulations on level zero or one seem to be a good compromise between 
complexity and efficiency. 

As shown by a variational calculation for high-energy scattering from a Gaussian potential the new repre- 
sentation offers new opportunities for approximations, seems to be applicable also at low energy and could be 
used for stochastic evaluation of real-time scattering similar as in Ref . |15j . At least it is an interesting addition 
to the standard body of knowledge in nonrelativistic scattering theory originating from field-theoretic methods. 
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